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Abstract 

Thermodynamic properties of the one-dimensional (ID) quantum 
well (QW) with miscellaneous permutations of the Dirichlet (D) and 
Neumann (N) boundary conditions (BCs) at its edges in the perpen¬ 
dicular to the surfaces electric field § are calculated. For the canonical 
ensemble, analytical expressions involving theta functions are found 
for the mean energy and heat capacity cy for the box with no applied 
voltage. Pronounced maximum accompanied by the adjacent min¬ 
imum of the specific heat dependence on the temperature T for the 
pure Neumann QW and their absence for other BCs are predicted and 
explained by the structure of the corresponding energy spectrum. Ap¬ 
plied field leads to the increase of the heat capacity and formation of 
the new or modification of the existing extrema what is qualitatively 
described by the influence of the associated electric potential. A re¬ 
markable feature of the Fermi grand canonical ensemble is, at any BC 
combination in zero fields, a salient maximum of cy observed on the 
T axis for one particle and its absence for any other number N of cor¬ 
puscles. Qualitative and quantitative explanation of this phenomenon 
employs the analysis of the chemical potential and its temperature de¬ 
pendence for different N. It is proved that critical temperature Ter of 
the Bose-Einstein (BE) condensation increases with the applied volt¬ 
age for any number of particles and for any BC permutation except 
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the ND case at small intensities S what is explained again by the 
modification by the held of the interrelated energies. It is shown that 
even for the temperatures smaller than the total dipole moment 
(P) may become negative for the quite moderate S. For either Fermi 


or BE system, the inhuence of the electric held on the heat capacity is 
shown to be suppressed with N growing. Different asymptotic cases 
of, e.g., the small and large temperatures and low and high voltages 
are derived analytically and explained physically. Parallels are drawn 
to the similar properties of the ID harmonic oscillator, and similarities 
and differences between them are discussed. 

1 Introduction 

The preceding paper [1] discovered, among other findings, the independence 
of the sign of the polarization on the boundary conditions (BCs) for 
the one-dimensional (ID) quantum well (QW) of the width L placed into 
the uniform electric field S that is directed perpendicular to its confining 
surfaces located at a: = ±L/2: the polarization Po(^) of the ground state for 
any permutation of the Dirichlet (D), 



( 1 ) 


and Neumann (N), 



( 2 ) 


edge requirements imposed on the wavefunction T(a;) is positive for all ap¬ 
plied voltages while its excited-state counterparts n > 1, for the 

small growing helds decrease from zero at ^ = 0 to the negative values, 
pass through the minimum and only after this start to increase crossing zero 
at the n- and BC-dependent intensity . Immediately, one wonders: for 
any kind of the particles, is it possible to observe the total statistically av¬ 
eraged polarization that is negative at the small electric forces? Analysis 
below answers this question together with the thermodynamic calculations 
of the corresponding energy E and heat capacity cy. Following the previous 
research [1], the QW with the particular distribution of the BCs will be de¬ 
noted by the two characters, where the first (second) one corresponds to the 
edge condition at the left (right) interface. Similar to the discussion of the 
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spectrum En{S') and polarizations Pn{^) [1], all energies will be measured, if 
not specified otherwise, in units of 7r^h^/(2mL^), which is a ground-state en¬ 
ergy of the DD QW, while the unit of the electric field will be (2emL^), 
and that of the polarization - eL, with m being the particle mass and e de¬ 
noting the absolute value of the electronic charge. In addition, heat capacity 
is expressed below in terms of Boltzmann constant ks- Discussion considers 
canonical as well as grand canonical ensembles. In this last case, the prop¬ 
erties are calculated both for fermions and bosons. Also, frequently we draw 
parallels with the ID harmonic oscillator (HO) with the potential (in regular 
units) [2] 


Vho{x) 


1 

= -mu 
2 


2 2 
X 


(3) 


whose energies upon application of the electric voltage, are 




hu 



2 mu‘^ 


(4) 


For this conhguration, the natural units that will be used below are: for the 
energy, tvu] for the length, xq = [h/(ma;)]^/^; for the electric held, hu/{exo); 
and for the polarization, exo- 


2 Canonical Ensemble 


This type of the statistical ensemble assumes that the system under consid¬ 
eration is in the thermal equilibrium with the much larger bath characterized 
by the thermodynamic temperature T. The fundamental quantity here is the 
partition function 

~ (5) 




/3-En 


where the summation runs over all possible quantum states, and the param¬ 
eter (3 is (in regular, unnormalized units) f3 = 1/(/cbT). The probability Wn 
of hnding particle in the state n depends on the temperature and the energy 
En as 

1 


— _ o — 


w„ = —e 


( 6 ) 


As a result, the mean value (X) 

1 


of any physical quantity X is calculated as 

(7) 


(X)can = ^ ^ 

-M ' ^ 


o — fiE-n 
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For the N particles in the system, this equation has to be multiplied by N. 
Applying these general results to the QW with the different BCs in the elec¬ 
tric field S’, one derives the mean values of the energy (E) and polarization 

(P) 


{P)can{f3,S) = 


— pEn 


EZpEne 

p-0Er, 

Z^n=0^ 
Z^n=0 
Z^n=0 


— 0En 


(8a) 

(8b) 


where in the left-hand side we have explicitly underlined that they are func¬ 
tions of the temperature T (through the parameter (3) and electric field 
[through the corresponding dependence of En{S) and Pn{S)]. Equivalently, 
Eq. (8a) can be written as: 


{E)can = -^InZ. 


(9) 


Heat capacity at the constant volume cy is a work that has to be done 
to change the temperature of the system by one degree and, as a result of 
this, it is calculated as a derivative of the total energy with respect to the 
temperature T: 

Cv = A(i5) = (10) 

where regular, unnormalized units have been used. Applying this generic 
definition to the canonical distribution from Eq. (8a), one gets fluctuation- 
dissipation theorem [2] 


(B>L) . 


( 11 ) 


where, for convenience of the notation, the subscript V has been dropped. 
Energies En and polarizations for the QW were calculated before [1] while 
for the HO they are: 


E^^ = n + ---S^ 
" 2 2 


pHO ^ _ 


dS 


= s. 


(12a) 

(12b) 


Note that, contrary to the hard-wall QW [1], for its HO counterpart the 
polarization is at any voltage a linear function of the field and is the same 
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for all levels. Accordingly, its mean value for the one particle is equal to S 
too while the energy becomes: 



can 





(13) 


As a result, the electric held does not affect the HO canonical heat capacity, 
which reads [2]: 

(14) 

One can derive limiting cases of these dependencies: 
for the small temperatures (/3 ^ cxo): 

i + + + + (15a) 

C = /Sqa-'* + 2a-'’+ 3a-W + ...) . (15b) 


for the large temperatures (/3 —?■ 0): 


\ / C 


^HO 


^ 1 1 1 
/3 ^ 12 ^ 720 

= 1-—/32 + — 
12 ^ 240 
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(16a) 

(16b) 


Before discussing the electric held inhuence on the thermodynamic prop¬ 
erties of the hard-wall QW, let us address hrst the voltage-free conhguration. 
Plugging in the well known expressions for the zero-held energies 


£""(0) = (n + 1)’, E^°(0) = (" + j) . "(O) = (17) 

into Eq. (9), one gets after some algebra: 


1 - 03 ( 0 , e-/5) d/3 ^ 

(18a) 


(18b) 


(18c) 
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Figure 1: (a) Heat capacity cy and (b) mean energy {E) as a function of the 
normalized temperature for the ^nonical ensemble and pure Dirichlet 
(dotted line), Neumann (solid curve) and ND (dashed line) QW at zero 
electric held. 












Here, 9i{z, q), i = 1,2, 3,4, are Theta functions [3,4], For small temperatures, 
/9 —)■ oo, these equations degenerate to 


(E^^) = 1 + 3e-3^ - 3e-®^ + + 36“^^ + ... (19a) 

{^^^)can = \ + - 2e"^^ + 8 e"®^ - 10e“®^ + ... (19b) 

/^JVTVX ^ g-/3 _ g-2/3 ^ g-3/3 ^ 3g-4/3 _ 4 ^- 5 /? 

\ / can 

+5e"®^ - 6 e“^^ + 3e“®^ + .... (19c) 

Utilizing transformation properties of the Theta functions [3] 

03 (O,e-^) = 03 ( 0 , 6 -^/^) (20a) 

O 2 ( 0 , e-P) = ^4 ( 0 , , ( 20 b) 


one derives the energies in the opposite limit of the high temperatures: 


En’^ 

/ can 


± 


/E^d) 

\ I c 


2/3 ^ 27rV2/3i/2 

1 27r2 


/S 


,-7r2//3 


/3 ^ 0 


-^2//3 


,5^0. 


(21a) 

(21b) 


The corresponding heat capacities cy are calculated by applying the right¬ 
most part of Eq. (10) to the above dependencies; in particular, one has 


for the “cold” QW, /3 —)■ cxo: 



= /32(gg-3/3_^gg-6/3^g4g-8/3^ 

,,,) 

(22a) 

= ^2(4g-2/3_gg-4/J^4gg-6/3_gQg-8/3 4.___) 

(22b) 

= /32 (e-^ - 2e-2^ + + 126' 

- 20e-®^ 


+ 30e"®^ -42e"^^-f 24e“®^+ ...) 


(22c) 

for the hot thermal bath: 



DD,NN _ 1 , 1 0I/2 „-vrV/3 

2 47rV2^ (3 ^ 

/3 ^ 0 

(23a) 

ND _ 1 47r2 2/^ 27r^ 2 /^ 

Ccan 2 /3 ^ ^ /32 ^ 

/3 ^ 0. 

(23b) 
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Figure 2: Mean energy {E) of the canonical ensemble in terms of the electric 
held S' and temperature for all permutations of the BCs. In each of the 
panels, the corresponding type of the edge requirements is denoted by the 
two characters. 







Statistically averaged energies and corresponding heat capacities are shown 
in Fig. 1. At the zero temperature, the total internal energy reduces to the 
ground-state energy and the heat capacity is zero, as it should be. As it 
follows from Eqs. (19) and (22), their avalanche growth from the T = 0 val¬ 
ues for the purely Neumann QW takes place at the smaller temperatures as 
compared to the mixed BCs, which, in turn, is followed by the quantities for 
the Dirichlet structure. This is explained by the growing difference between 
the two lowest energies for just consecutively mentioned BC conhgurations: 
the quantity 

= En+i{S) - En{S) (24) 

at the zero held is the smallest (largest) for the Neumann (Dirichlet) struc¬ 
ture: 


A(y^(0) =2n + l, A)y^(O) = 2n + 2, A^^(0) = 2n + 3. (25) 

A remarkable feature of the heat capacity dependence is its nonmonotonic 
behaviour for the Neumann QW: at = 0.4342 {Pmax = 2.3031) it reaches 
a pronounced maximum ~ 0.4455 that is followed by the minimum of 

^min = 0.3818 located at = 0.9420 (/9mm = 1.0616). If the maximum is 

observed quite exactly by keeping only the hrst term in the parentheses of 
the right-hand side of Eq. (22c), the emergence and precision of the location 
and magnitude of the second extremum are described better by keeping more 
terms in the same expansion. Physically, this nonmonotonicity of the heat 
capacity is attributed to the structure of the energy spectrum, see Eqs. (17) 
and (25); namely, very small temperature promotes the particle mainly to the 
hrst excited level that is only one unit above the ground state, A^^(O) = 1, 
with the contribution of the other levels being negligibly small due to the 
almost vanishing exponents in Eq. (8a) or, equivalently, in Eqs. (19c) and 
(22c); as a result, the heat capacity grows rapidly. For the larger tempera¬ 
tures, the occupations of the higher lying levels become essential; however, 
the transitions to them are more difficult since the difference between, e.g., 
second and hrst excited states A(^^(0) = 3 is three times larger than that 
between the latter and the ground level. Accordingly, the same speed of 
the heat capacity change can not be sustained what results in the observed 
maximum. For the other BCs, the ratio Ai(0)/Ao(0) is smaller than for 
the Neumann QW, as it follows from Eq. (25): A(^^(0)/A^^(0) = 2 and 
Af^(0)/A^^(0) = 5/3; as a result, for them no extrema are observed on 
the cy — T dependence at /9“^ < 1. Mathematically, the drop of the NN 
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specific heat is caused by the interplay between the counterbalancing terms 
and e~^ in Eq. (22c) as the temperature grows. Keeping only the hrst 
exponent in the parentheses of the right-hand side of this equation produces 
= 2 while the same procedure applied to the other BCs, see Eqs. (22a) 
and (22b), results in f3max^'’ =2/3 and (3max^'‘ = 1, which are, respectively, 
three and two times smaller and lie beyond the range of the validity of these 
expansions. Accordingly, for the latter two conhgurations, it is essential to 
keep other items in the corresponding series in order for them to be correct 
at the decreasing /3, and these extra exponents eliminate the resonance of the 
hrst-term approximation while for the Neumann QW the (negative) second 
component simply improves the previous result. Note that the HO leading 
term of the capacity expansion from Eq. (15b) also results in I3max^'’ = 2; how¬ 
ever, the subsequent (all positive) items in the series wipe out the extremum. 
Very broad and gentle asymmetric maximum is observed at > 2.5 for the 
Dirichlet QW while for the mixed BC the heat capacity is a monotonically 
increasing function of the temperature, which, at quite large T, rapidly ap¬ 
proaches the asymptotic value of one half. On the contrary, the heat capacity 
of the symmetric QWs reaches the same limit much slower, as Eq. (23a) as¬ 
serts and panel (a) of Fig. 1 exemplihes. Note that the HO internal energy 
for the high temperatures is twice of that for the hard-wall QW: ksT and 
^ksT in regular units, respectively. From point of view of classical equilib¬ 
rium statistics that is applicable for T —)■ oo, this difference is explained by 
the fact that in the former case the kinetic and potential parts of the motion 
make equal contributions of ^ksT to the total energy [2] while for the latter 
system it is the kinetic energy only that determines {E) as the QW potential 
is zero. As a direct consequence of this, the QW heat capacities in the same 
limit are one half of their HO counterpart. 

Applied electric held modihes the energy spectrum what, in turn, affects 
the thermodynamic properties of the wells. It was shown that the voltage 
increases the difference Aq{S‘) between the ground and hrst excited levels 
for any permutation of the BCs (the only exception is the ND case at the 
small helds, see equations (50) in [1]); accordingly, the larger temperature 
is needed to push out the electron from its lowest state. This is rehected in 
Figs. 2 and 3 where the energy {E{(3, S‘))can and heat capacity {cv{{3, <f)}can, 
respectively, are shown. It is seen that the range where the mean energy 
does not change appreciably from the ground-state value gets wider for the 
stronger intensities S. The same is true for the heat capacity where the 
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Figure 3: The same as in Fig. 2 but for the heat capacity cy. 
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plateau with its almost zero value grows with the held. The increasing voltage 
wipes out the NN minimum of the heat capacity simultaneously moving the 
maximum to the higher temperatures and increasing its magnitude. For 
each of the mixed BCs, it also creates a maximum that was absent at ^ = 0. 
Mentioned above DD extremum of the heat capacity gets narrower and its 
peak increases with the held growing. Recalling again the language of the 
classical statistical mechanics [2], one qualitatively explains the larger heat 
capacities at the nonzero helds by the contribution of the electric potential; 
namely, the thermally averaged value of the potential energy {—Sx) is: 

(-.Ti) = i b - 1/3.? coth i/3#V (26) 

This classical expression is applicable to our quantum system for the large 
temperatures only: 

{-■S'x) « - ..., /3 ^ 0. (27) 

Then, the potential contribution to the heat capacity reads: 

_ ^(f).r)4 +..., (28) 

Note that, contrary to the HO, the kinetic and potential contributions to 
the heat capacity in this case, generally, are not equal to each other. Let us 
also mention once again that the electric field does not affect at all the HO 
heat capacity, see Eq. (14), since it simply shifts all the levels by the same 
amount, according to Eq. (12a). 

Fig. 4 depicts statistically averaged polarizations {P)can hi terms of S’ 
and Growing temperature leads to the decrease of {P)can for all electric 
helds; however, thermal energy is not strong enough to make the total dipole 
moment negative: for any BC the polarization stays positive. To understand 
the statistical properties better, it is instructive to consider the case of the low 
temperatures. For the small voltages, as a hrst approximation, we also accept 
undisturbed by the held energies from Eq. (17). Then, one has following 
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Figure 4: The same as in Fig. 2 but for the polarization (P). 
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dependencies: 


^pDD) 

= 

Po + (Pi - Po) {e-^^ 

- e-6^) 



+ 

(P 2 -Po)e-®^ + ..., 

/3 —)■ oo 

(29a) 

^pND) 

= 

Po + (Pi - Po) (e-2^ 

1 

1 



+ 

(P2 - 2Po + Pi) 

-|- ..., /3 —)■ CX) 

(29b) 

(^pNN^ 

= 

Po + (Pi - Po) (e-^ - 




+ 

(P2-Pi)e-"^ + ..., 

/5 —)■ cxo. 

(29c) 


These equations were derived under the assumption of <C 1, but the general 
property stating that the hrst-order temperature correction is determined by 
Pi — Pq, holds for any electric intensities. For the small helds, this difference 
is negative [1] what naturally explains the decrease of the total polarization 
with the temperature growing. In the opposite limit of the high voltages, 
the polarizations of the QW with the uniform BCs tend to the same level- 
independent value of one-half [1] what requires larger temperatures in order 
to see the deviation of (P) from its T = 0 value. This is exemplihed in Fig. 4 
where the temperature-independent plateau at T = 0 widens with the held 
growing. For the mixed edge requirements, this limiting quantity is supple¬ 
mented by the term that is proportional to ±(n -|- 1/2)“^ with its sign being 
determined by the orientation of the BCs [1]; so, for the DN case it is actu¬ 
ally possible to observe the increase of the polarization with the temperature 
growing from zero. This feature is not shown in the corresponding panel of 
the hgure since it takes place beyond the hgure range 0 < (^ < 50. 


3 Grand Canonical Ensemble 

Grand canonical distribution is used for the description of the quantum sys¬ 
tem that, in addition to the thermal balance with the external reservoir, is 
also in the chemical equilibrium with it. Accordingly, the structure can ex¬ 
change the energy as well as particles with the heat bath. So, the number 
of the quantum corpuscles N in it can be changed. The fundamental role 
in this case is played by the chemical potential p, which is dehned from the 
condition 

^ ^ X] ± 1 ’ 
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where the upper sign corresponds to the Fermi-Dirac (FD) distribution while 
the lower one describes Bose-Einstein (BE) particles. Physically, the differ¬ 
ence between these two statistics is in the fact that each quantum level can 
not be occupied by more than one fermion (Pauli exclusion principle) while 
the arbitrary number of bosons can coexist in the same state. The distri¬ 
bution function now depends not only on the energies but also on the 
number of the particles in the system iV; namely, for the physical quantity 
X its grand canonical average value (X)gc is 

(^)5c = e(En-/.)/3 ± 1 • 


Applying the distribution from Eq. (31) for the calculation of the heat ca¬ 
pacity, Eq. (10), one hnds that its grand canonical value Cgc is 


Cgc 




E„ - /I - /3§ 




(32) 


where the chemical potential, which, in the case of fermions, is also frequently 
called the Fermi level, is calculated, as stated above, from Eq. (30). Phys¬ 
ically, the value of p corresponds to the energy that is needed for changing 
by one the number of the particles in the system: 





(33) 


For calculating its partial derivative with respect to the temperature, one 
should consider Eq. (30) as a condition of zeroing of the implicit function 
F(/i, /9, N) of the chemical potential in terms of the variables (3 and N\ 


F{^^{(3,N),(3,N) = Q. 


The rule of differentiating implicit functions states [9]: 

dF/d^ 

^ “ ~dFld[i' 

As a result, one hnds: 



E 


_ ^ g(£^7T,—/7,)/3 

g(Sn —m)/3 

-1- 

[e^Er,-^)P±lY 


(34) 


(35) 


( 36 ) 


16 










The boson statistics is used for the particles with the integer spin such 
as photons or Cooper pairs in superconductors while the FD distribution is 
applied for the system of the constituents with the half integer spins; for 
example, the electron with its spin of 1/2 has, for the same energy, two 
projections of its spin equal to ±1/2. However, in our discussion below we 
will neglect this fact and will assume that the number of the fermions for 
each energy is not larger than one. 

Fig. 5 depicts the FD chemical potential for the pure Dirichlet QW as a 
function of the electric held and temperature for several numbers N. Quali¬ 
tatively, the same features are characteristic for other BC permutations too. 
There are several distinct regions of the Fermi energy dependence on the 
temperature. From its Ej^_i value at T = 0 it rapidly grows as 




En-i — 



1 ± \/l - 4e-^Jv-i/3 
2 


—)■ cxo. 


(37) 


until it reaches and stays exactly at the value of 

hAf = 2 {Em-1 ± Em) , ^ 1, 


(38) 


which is due to the interaction of the two corresponding levels that, at the 
zero temperature, were the highest occupied and lowest unoccupied states. 
The width of this T-independent plateau is determined by the number of the 
particles N and electric held S. For the still higher temperatures, the chem¬ 
ical potential for = 1 decreases while for the larger number of particles, 
> 2, it grows with reaches maximum and only after that decreases, 
passes zero at and continues to decline into the negative part of the 
spectrum. For the nonpositive chemical potentials, \im < 0, Eq. (30) can be 
cast into the form 


OO OO 

^- 0 {m+l)En _ jy 

m=0 n=0 


/i < 0. 


(39) 


For the zero held, ^ = 0, the energy spectrum from Eqs. (12a) and (17) 
simplihes this equation as follows: 
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Figure 5: Chemical potential for the FD distribution of the pure Dirichlet 
QW as a function of the electric held S’ and temperature for different 
number of fermions N that is depicted in each panel. Note different /i scale 
for each plot. Viewing perspectives for the upper two panels are different 
from those for the other subplots. 
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for the HO: 


OO 


1 

2 


Efri)' 


g/i^(m+l) 


m=0 


sinh i/3(m + 1) 


= N, /i < 0 


and, for example, for the mixed BCs: 

OO 

^(^l)mg/./3(m+l)^^ (0, =N, /i < 0. 


m=0 


(40a) 


(40b) 


Putting here the chemical potential equal to zero, p = 0, leads to the cal¬ 
culation of In known to us literature [3-8], there are no analytical 

expressions for these inhnite series. However, for the very small /9, the m = 0 
terms in the above equations make the most signihcant contributions pro¬ 
ducing the following dependencies: 
for the HO: 

= (41a) 

and, for the hard-wall QW with the arbitrary BCs: 


Superscripts I and J in Eq. (41b) stand for any of the values of D and/or 
N. Fig. 5 manifests that, for the larger N, these asymptotics are achieved 
at the higher temperatures. As a result, the grand canonical mean energy 
{E)gc reads in the same limit: 

= 0^0 (42a) 

(•E'''U/3)U.„ = ~. (3^0, (42b) 

what, by means of Eq. (10), immediately leads to the associated heat capac¬ 
ities Cgc'. 

4°(/5)C„ = 3V, P^O (43a) 

4iW\i!.o = Y’ 

A comparison of these remarkable results with Eqs. (16), (21) and (23) con- 
hrms the general property, which states that for the large temperatures there 
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is no difference between canonical and grand canonical distributions [2]. How¬ 
ever, for the small T these two statistics produce very different features. 
Fig. 6 shows the FD heat capacity of the pure Dirichlet QW in terms of 
the temperature and electric field for the different N corresponding to their 
counterparts from Fig. 5. It is seen that, for the larger number of the par¬ 
ticles, the asymptotics from Eq. (43b) is achieved at the higher T. At the 
zero field, a prominent characteristic of the heat capacity dependence for the 
one particle (top left panel of Fig. 6) is a salient maximum Cmax = 0.882 
observed at = 0.633, i.e., at the right edge of the plateau from Eq. (38). 
Accordingly, we attribute this extremum to the different behaviour of the 
chemical potential for A^ = 1 and N > 2] namely, as it was mentioned during 
discussion of Fig. 5, for one particle the Fermi energy decreases after the flat 
part from Eq. (38) while for any other number N it grows with T. Thus, their 
contributions to the heat capacity from Eq. (32) are opposite to each other 
what results in the resonance that is observed for the one particle only. Even 
though the shape of this maximum is quite similar to its NN counterpart 
for the canonical ensemble, see Sec. 2, its physical explanation is completely 
different. First, we point out that the very similar extrema are calculated 
also for the ND (with Cmax = 0.879 and = 0.418) and pure Neumann 
{cmax = 0.878 and P^ax = 0.208) QWs too. The fact that the three Cmax are 
almost the same and the ratios of the three temperatures T^ax are practically 
equal to those of Ao(0) from Eq. (25), undoubtedly proves that the origin of 
this effect is the BC independent one and that the interplay between the two 
lowest states plays a dominant role in it. To understand these resonances, let 
us recall that, for the very small temperatures, the properties of the FD well 
are determined only by the highest occupied level and its interaction with the 
nearest (empty at T = 0) above lying state, what is reflected in the extremely 
rapid approach by the chemical potential to the energy from Eq. (38) that 
is located exactly in the middle between them. For N > 2, a contribution 
from the lower lying members in this regime is negligibly small and can be 
safely neglected, while for the one-electron QW this addition is absent by 
dehnition. Further growth of the temperature increases thermal energy but 
it is still too “weak” to compel the corpuscles, which at T = 0 lied below the 
Fermi energy, to contribute to the heat capacity. Only at the right edge of 
the plateau, the thermodynamic quantum ksT becomes strong enough and 
forces other particles to donate to /tat and cn- Therefore, for A^ > 2 the heat 
capacity is a quite smoothly varying function of the temperature. However, 
for N = 1 there are no such additional donors that aid to support the con- 
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tinuous growth of the heat capacity, which can not be sustained by the one 
particle only. As a result, the specific heat reaches maximum and drops. 

This qualitative physical reasoning can be corroborated by the simple 
quantitative mathematical analysis. Fermi level from Eq. (37) defines the 
corresponding mean energy and heat capacity as: 


{E) = J2En + lEN-i + EMe-^^-^^ 13^ 00 (44a) 

cy = /5 —)■ cxo. (44b) 


On the other hand, for the chemical potential from Eq. (38) these quantities 
for one fermion, = 1, become: 


(E) 


Cv 


En 


+ 


El 


1 + (>-AoPI2 ]_ gAo/3/2 ’ 


/5>1 




Er 


eAo/3/2 


(1 + e^o/5/2) 
g-Ao/3/2 


(1 + 


/5>1, 


(45a) 


(45b) 


where we 
Neumann 
pie form: 


take into consideration only the two lowest states. For the pure 
QW without the field, this last expression takes an especially sim- 


NN\ 

1 1^=0 
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(46) 


This function has a pronounced maximum of Cmax = 0.878 at Pmax = 4.799. A 
perfect coincidence with the provided above exact results justifies a validity of 
the two-level approximation and proves that the electron transitions between 
them determine the specihc heat resonance. It is also very instructive to 
contrast Eq. (46) with its canonical counterpart from Eq. (22c) in sec. 2. 
The comparison shows that the magnitude of the grand canonical Neumann 
extremum is almost two times larger and it is achieved at more than two 
times lower temperature. 

Applied held S’ smooths out and widens this maximum simultaneously 
increasing the heat capacity. Similar to the canonical ensemble, this growth 
is explained by the contribution of the electrostatic potential. However, the 
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Figure 6: The same as in Fig. 5 but for the heat capacity c^. Note different 
c and scales for different panels. 
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electric influence is drastically decreased by the growing number of the par¬ 
ticles in the QW; for example, right-bottom panels exhibit the almost full 
independence of the Fermi energy and heat capacity on the intensity S al¬ 
ready for N = 10. This is explained by the properties of the energy spectrum 
in the electric field when the higher lying states (which, in the case of the 
FD distribution, determine the features of the system) are less affected by 
the applied voltage [1]. 

Fig. 7 demonstrates zero-temperature FD polarization for all possible 
BCs and several numbers iV. As the well accommodates more fermions, the 
total polarization becomes smoother function of the electric held. Figure 
reveals that, independently of the edge demands, the magnitude of (P) at 
N > 5 grows linearly with the voltage and the slope of this almost straight 
line diminishes with N. For any number of fermions, the total polarization 
remains positive at the arbitrary voltage. Nonzero temperature leads to the 
dependencies that qualitatively are similar to the canonical patterns. Fig. 4, 
and, because of this, the corresponding polarizations are not shown here. 

Next, let us discuss bosonic structures. Remarkable experimental obser¬ 
vations of the BE condensation in the vapors of rubidium [10] and sodium 
[11] spurred an avalanche of the research on the subject predicted almost 
ninety years ago [12], see, e. g., reviews [13-16]. Theoretically, the main 
effort was devoted to the calculation of the properties of the BE systems 
in the 3D isotropic or anisotropic harmonic traps [13,17-20] and their ex¬ 
istence/nonexistence in lower dimensions [13,17,18,20,21]. However, other 
forms of the confining potentials [21,22], including the 3D box with the 
periodic [23-25] or uniform [22,26-29] BCs, were also discussed with the 
comparative analysis of their influence of the properties of the trap [30-33]. 
From this point of view, an inclusion of the electric voltage and different BCs 
presents a generalization of the previous analysis. Moreover, overwhelming 
majority of the research concentrated on the analysis of the BE systems in 
the thermodynamic limit when the number of the particles and the volume 
containing them tend to infinity while the the density is kept constant. In 
this approximation, the inhnite series above in this section can be safely re¬ 
placed by the integrals [13,14,16]. Considering N changing from one to the 
large values might help to understand the formation of the BE processes with 
the the number of the particles growing. First, we state that Eqs. (41) - (43) 
stay valid for the BE statistics too since they were obtained as a result of 
retaining the first term in the series from Eq. (40). In the opposite limit of 
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Figure 7: Total polarization (P) as a function of the electric field S at zero 
temperature, T = 0, for different number of fermions N that is depicted next 
to the corresponding curve. 
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the very small temperatures, it is elementary to derive: 


Pjv = ^0 - ^ , /3 ^ oo, (47a) 

what leads to the mean energy and heat capacity: 

{E)n = NEo + (3 ^ oo, (47b) 

cn = ^ EiAo(3‘^e~^°^, (3 -i- oo. (47c) 


An important characteristic of the BE system is its critical temperature 
Ter- It corresponds to the situation when the chemical potential is equal to 
the energy of the lowest level, fi = Eq, and the number of the particles in 
this state Nq is zero what leads to the implicit mathematical equation for 
hnding (3cr 


E 

n=l 


1 

fp{En—Eo)lii 

cr 


1 


N. 


(48) 


Physically, it is the largest temperature at which the BE condensation still 
can be observed, and at the lower T the fraction N^/N of the particles in the 
ground state will increase until at T = 0 it becomes unity: 


A'o/A'It=o = 1- («) 

Fig. 8(a) shows dependencies of the critical parameter on the applied 
voltage for all possible BCs and several numbers N. In accordance with the 
previous results [13], the temperature Ter increases with N. In the absence 
of the helds, the lowest (highest) temperature is observed for the pure Neu¬ 
mann (Dirichlet) QW what is a reflection of the corresponding spectrum from 
Eq. (17) and the energy difference between the affiliated states, see Eq. (25). 
Electric field leads to the modification of the mutual location of the levels 
on the energy axis; in particular, at the small voltages, the two lowest states 
move closer to each other for the ND geometry while the difference Ei — Eq 
grows with the field for all S and any other BC configuration [1]. As a result, 
the critical temperature for the former edge requirement decreases with the 
growing from zero field, passes through minimum and then unrestrictedly 
grows with the electric intensity while for all other BCs it is a continuously 
increasing function of S. At the high voltages, the energy spectrum is de¬ 
termined mainly by the condition at the right wall [1] what leads to almost 
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Figure 8: (a) Critical temperatures and (b) corresponding to them po¬ 
larizations {Per) from Eq. (50) as fictions of the applied electric held S 
for different number N of bosons that are depicted near the corresponding 
curves. Solid (dotted) lines are for the pure Dirichlet (Neumann) BC while 
their dash (dash-dotted) counterparts denote DN (ND) geometry. In panel 
(b), thin horizontal line is zero polarization. 














the same critical temperature for, e. g., the NN and DN wells. Note that 
in this regime, contrary to the zero helds, the Dirichlet requirement is more 
favorable to the formation of the BE condensate as compared to the Neu¬ 
mann interface. This is explained by the larger level separation for the latter 
geometry at the high voltages [1]. 

As the ground-state polarization is positive for all fields and BCs [1], one 
can expect that at the onset of the BE condensation the total statistically 
averaged dipole moment will, at the small voltages, be negative. This is 
exactly what is observed in panel (b) of Fig. 8 that shows the polarizations 
{Per) corresponding to the critical temperatures from panel (a). They 
were calculated from equation 


p 

{Per) = ^ ^{E^-EQ)fl,r _ 1 ’ 
n=l 

and, as stated above, the critical temperature was found from Eq. (48). 
The characteristic features of the critical polarizations basically follow the 
properties of the first excited state: from zero they decrease with the growth 
of the held, reach minimum after which they increase. However, for the large 
voltages many upper lying states are occupied and contribute to the total 
dipole moment. As a result, the high-held {Per) for one boson is considerably 
smaller than Pi. The absolute value of the negative polarization at the 
extremum grows with the number of the particles with the largest one, at 
the hxed iV, being observed for the pure Neumann QW followed by its ND 
counterpart what is a replica of the similar behaviour for the hrst excited 
level [1]. For the temperatures below Ter, the nonzero occupation of the 
ground state contributes a positive term to the polarization what leads to 
the gradual disappearance of the negative region of the total dipole moment 
(F) with the decreasing temperature until at T = 0 it becomes NPq, which 
is positive for any BC and arbitrary helds [1]. 

As a hnal example. Fig. 9 exhibits evolution of the heat capacity and 
chemical potential with the varying electric held and temperature for the 
pure Neumann QW. It is seen that for the small number of bosons, say, 
= 3 in panel (a), the applied voltage leads, at quite warm sample, to the 
increase of cy while at the small T, the width of the temperature-independent 
zero-capacity plateau increases with the held. These features were discussed 
before for the canonical ensemble. Increasing the number of the particles 
in the well leads to the suppression of the voltage dependence, as a tran- 
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Figure 9: Bosonic heat capacity cy of the pure Neumann QW as a function of 
the applied electric held S and temperature for (a) iV = 3, (b) iV = 1000 
and (c) N = 10®. Note that in the last two cases the temperature is scaled in 
units of the critical temperature while the corresponding 2 ; axes measure 
specihc heat per particle c/N . Panel (d) shows the chemical potential fi for 
N = 10® with the corresponding heat capacity depicted in part (c). 
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sition from panel (a) to (b) with N = 1000 and (c) for N = 10^ vividly 
demonstrates. No any noticeable field dependence is seen there in the range 
0 < (a < 50. It is well known that for some potentials, such as, e.g., the 3D 
isotropic harmonic trap [14,16,18,19,34,35], the heat capacity has a cusp¬ 
like peculiarity as it passes through the critical temperature while for the ID 
quadratic potential it is a smooth function of T [16,18,35]. Fig. 9 exempli¬ 
fies that no any peculiarity is observed for the ID hard-wall potential with 
Neumann surfaces and arbitrary applied electric helds. Our calculations con¬ 
firm that the same is true for any other BCs. Finally, panel (d) shows that 
the chemical potential p is a monotonically decreasing function of both the 
electric held (p and temperature T. It is seen that the growing temperature 
diminishes the voltage inhuence on the chemical potential. 


4 Concluding remarks 

Rigorous mathematical treatment of the QW with miscellaneous BCs under 
the applied voltage revealed a strong inhuence of the interplay between them 
on the thermodynamic properties of the structure. In particular, without the 
held the diherences of the energy spectrum lead, for the canonical ensemble, 
to the conspicuous maximum followed by the minimum of the heat capacity 
cy on the temperature axis for the NN quantum box while for the other 
edge requirements no such adjacent extrema are observed. Modihcation of 
the specihc heat and statistically averaged polarization in the held is quali¬ 
tatively explained by the inhuence of the associated electrostatic potential. 
Numerical calculations, which predicted, for the hat potential with the arbi¬ 
trary BC, a salient maximum of cy as a function of T for one fermion and 
its absence for the larger N, were corroborated by the two-level model that 
allows simple analytical treatment with its predictions perfectly coinciding 
with the exact results. From this, a clear physical explanation of this phe¬ 
nomenon follows that is based on the analysis of the associated Fermi energy. 
It is predicted that the applied held, in general, favors the formation of the 
BE condensate, and the diherences and similarities of this process for the dif¬ 
ferent BCs are discussed. The thermally averaged dipole moment is shown 
to take the negative values in some ranges of the helds and temperatures. 
It is also argued that for the larger number of either fermions or bosons in 
the QW, the inhuence of the electric held on the thermodynamic properties 
diminishes. 


Dirichlet and Neumann conditions are the limiting cases of the so called 
Robin BC [36] 

> ( 51 ) 

s 

where n is an inward unit normal to the surface, and the parameter A has a 
dimension of length and is called the extrapolation distance. Its variation al¬ 
lows a continuous transformation from the Dirichlet (A = 0) to the Neumann 
(A = oo) situation. Without the held, especially intriguing are the properties 
of the QW at the small negative Robin lengths, A —0, when, in addition 
to the positive spectrum, two almost degenerate odd and even levels with the 
energies E ~ —l/(7rA)^ are created [37]. Analysis of the Robin QW in the 
electric held might present an interesting extension of the present research. 
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